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Abstract. In this series of lectures, I shall begin with the current investigations on 
phenomenology of hadron dynamics to demonstrate the importance of solving hadronic 
bound states within the framework of light-front (LF) QCD. Then, I will describe the 
basic procedure how to formulate the canonical theory of LFQCD, including light- 
front quantization of QCD, light-front gauge singularity, and light-front two-component 
formalism. I will also present a complete one-loop QCD calculation in terms of the 
light-front time-ordering perturbation theory, in comparison with the usual covariant 
perturbative QCD calculation. Following thereby I will discuss the development of 
heavy-quark effective theory and the manifestation of heavy quark symmetry on the 
light-front. Finally, by applying recently developed similarity renormalization group 
approach to light-front heavy quark effective theory, I will show a rigorous derivation of 
quark confinement interaction from LFQCD and its application to solve heavy hadron 
bound states. 



1 Hadronic Phenomenology in the LF Formulation 
1.1 An Overview 

Simply speaking, the main task in the investigation of hadronic physics is how to 
provide a QCD description of hadronic structure. More specifically, how can we 
compute directly from QCD the fruitful hadronic properties, such as the hadronic 
structure functions in lepton-nucleon deep inelastic scatterings, the partonic frag- 
mentation functions in high energy hadron-hadron or e^e~ cohisions, and many 
hadronic form factors in various hadronic decay processes. However, although 
QCD has been accepted as a fundamental theory of the strong interaction that 
governs the underlying dynamics of hadronic constituents, a complete QCD de- 
scription to hadronic structure is still lacking. In this series of lectures, I will 
attempt to show you that the light-front formulation of field theory may provide 
a natural and systematic QCD description to all the processes mentioned above 
[Zhang (1994)]. 

Historically, light-front dynamics played a very important role in every step 
of the development of the strong interaction theory. The most important appli- 
cation of light-front dynamics to hadronic physics is perhaps the parton phenom- 
ena in the lepton-nucleon deep inelastic scatterings (DIS). As it is well-known, 
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DIS probes hadronic dynamics near the light-cone. Physically, the DIS phenom- 
ena can be understood in terms of Feynman's parton picture [Feynman (1972)]. 
While, only the light-front formulation of field theory can provide a natural 
quantum field theory description of parton dynamics. For examples, the leading 
contributions of the unpolarized structure function ^2(2;, Q^) and the polarized 
structure function gi{x, Q^) in DIS are simply written in terms of hadronic ma- 
trix elements on the light-front surface ^~ = 0: 

^^^^^^ = ^ / rfr/e-"'-(ps|^;(r)Q'V'+(0) - h.c.\ps), (1) 

9iix,Q^) = ^ J rf»7e-'''^(ps|vUr)75QV+(0) + h.c\ps). (2) 

Here, = 57°7^V-'(C) is the light-front quark field operator, 77 = ^p'^^~, 

^~ = — C^, Q the quark charge operator, and \ps) the hadronic states. It can 
be shown that in the light-front field theory, Eqs.(l) and (2) are proportional to 
the momentum and helicity distributions of partons (quarks and gluons) inside 
hadrons respectively. Other structure functions {Fl and (72) also have a similar 
but a bit complicated expressions. Nevertheless, it is obvious that if we knew the 
hadronic bound states \ps) from QCD on the light-front, we could completely 
understand the QCD dynamics of DIS. 

Another measurement of hadronic structure in terms of light-front hadronic 
matrix elements is the parton fragmentation functions in hadron-hadron and 
other collisions. During high-energy collisions, many hard partons are produced 
and then arc hadronized. Hadronization processes can be characterized by the 
so-called fragmentation functions which is also introduced initially by Feyn- 
man [Feynman (1972)]. Physically, quark and gluon fragnicintation functions are 
probabilities of finding hadrons in a hard parton produced in collisions. These 
fragmentation functions can be defined as matrix elements of quark and gluon 
operators at light-front separations. For examples, the unrenormalized quark 
fragmentation function is given by [Collins and Soper (1982)] 

fA/,iz) = dx-e-'^^^-/^Tr{0\ij+{0)\ps)(j>s\i,l{x-)\0), (3) 

and the gluon fragmentation function is defined as 

■^^/^^^^ = 36^+ / dx-e-"^'^-/^0\F+^mps){ps\F+{x-m. (4) 

where Tr traces the color and Dirac components of quarks. Again, if we knew 
hadronic bound states from QCD on the light-front, wo could directly study the 
QCD dynamics represented by these fragmentation functions. 

In recent years, light-front formulation has also been widely used in the phe- 
nomenological study of hadronic form factors involving in various hadron elastic 
scatterings and decay processes, by the use of the so-called relativistic quark 
model or light-front quark model [Terent'ev et al. (1976)]. Simply speaking. 
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light- front quark model is based on truncated Fock space expansion of light-front 
bound states (upon only the valence quark states) and then phenomenologically 
determines the valence Fock states' amplitude (the wavefunction). Unlike the 
study of the structure functions and the fragmentation functions where the use 
of light-front description can make the physical picture manifestation, the inter- 
esting feature of using light-front description to hadronic decay processes is that 
the simple boost operations and the transparent relativistic properties contain- 
ing in light-front bound states may allow one to describe hadronic form factors 
for entire kinematic range of momentum transfer for these space-like processes. 
This is quite different from descriptions of other hadronic quark models, such as 
the nonrelativistic constituent quark model and the beg model, which are nor- 
mally believed to be applicable only for the processes involving small momentum 
transfer. Very recently, applic;ations of light-front quark model have also been 
extended to the description of various heavy meson decay processes, although 
most of the investigations are limited to the calculations of form factors at zero 
momentum transfer, due to the limitation of using the light-front quark model 
for time-like processes. Extending the light-front quark model incorporated with 
higher Fock space contribution (a more realistic light-front bound state descrip- 
tion) may make the description of hadronic decay form factors become possi- 
ble for the entire kinematic range of momentum transfer. Nevertheless, all the 
hadronic form factors are extracted from some hadronic matrix elements, such 
as, 

{H'ip')\r\H{p)), (5) 

where F is a transition operator in the corresponding process. Again, if we knew 
the associated hadronic bound states that solved from QCD on the light-front, 
wc would have a true QCD description of hadronic decay processes. 

The above analysis indicates that once we know how to solve the hadronic 
bound states from QCD, especially for these defined on a surface of light-front, 
we can directly calculate various hadronic matrix elements involved in many 
hadronic processes. Then a true QCD description of hadronic physics may be 
realized. This series of lectures is devoted to the light-front formulation of QCD 
dynamics and the attempt of solving hadronic bound states, especially the heavy 
hadron bound states, directly from such a formulation. In the first Icctiire, I will 
mainly discuss the general structure of hadronic bound states on the light-front. 

1.2 General Structure of Light- Front Bound States 

In the standard language of field theory, relativistic bound states and resonances 
are identified by the occurrence of poles in Green functions. Although the in- 
formation extracted from this approach provides a good definition of physical 
particles, the ordinary wave function structure of bound states in the usual 
quantum mechanics language is lacked. As a result, wave function amplitudes 
extracted from Green functions may not be universally valid in the calculations 
of various hadronic matrix elements that measured in experiments. In order to 
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understand hadronic structure in terms of hadronic bound state wavefunctions 
(which is the most transparent picture in quantum theory), the expUcit form of 
hadronic bound states on some fixed time surface is wanted. 

However, solving bound states in field theory as an cigcnstatc problem has 
not been well established. One may define the bound states as eigenstates of 

and determine these states by solving the eigcncquation of P°. But P° is 
a square root function of the momentum and mass operators which does not 
give us a clear picture of the Schrodinger's eigenstate equation in quantum me- 
chanics. The widely used framework of finding relativistic bound states is the 
Bethe-Salpctcr equation. However, Bcthc-Salpctcr equation itself involve many 
unsolved problems, such as the physical interpretation of the Bethe-Salpeter am- 
plitudes, and the numerical difficulty in solving the Bethe-Salpeter equation in 
space-time space;, etc. Some approximations, such as instant-time approxima- 
tion, may simplify the Bethe-Salpeter equation. But with such approximation, 
the main properties of relativistic dynamics, namely the boost dynamics, will be 
lacking. In other words, the results may be no longer relativistic. 

Also, in principle, a relativistic bound state can always be written as an 
operator function of the particle creation operators acting on the vacuum of the 
theory. However, for many theories that we are interested in, especially for QCD, 
the vacuum is very complicated. With a complicated vacuum, formally writing 
down a relativistic bound state as a series of Fock space expansion also becomes 
very difficult. 

However, these subtle problems may be removed when we look at the bound 
states on the light-front. 

i) . Light-Front Vacuum. In the equal-time framework, the vacuum of QCD is 
crucial for a realization of chiral symmetry breaking and color confinement. It 
is also a starting point in the construction of hadronic bound states. However, 
the understanding of the true QCD vacuum is still very limited, although a 
lot of informative work has been carried out in the past two decades based on 
the instanton phenomena [Hooft (1976)] and the QCD sum rule [Shifman et al. 
(1979)]. 

In the light-front coordinates, a particle's momentum is divided into the 
longitudinal component and the transverse components. For a physical (on-mass- 
shell) particle, its longitudinal momentum, = kP -\- , cannot be negative 
since the energy of a physical state always dominates its momentum. As a result, 
the light-front vacuum for any interacting field theory can only be occupied by 
the particles with zero-longitudinal momentum, namely 

\vac)LF = /(4+=o)|0) > (6) 

so that P^\vac) LF = 0, where = '^ik^ . At this point, the light-front 
vacuum is still not simple. In the past several years, to obtain a nontrivial light- 
front vacuum, many tried to solve the so-called zero-mode (the particles with 
k+ = 0) problem [Burhkart (1996)]. 

To construct hadronic bound states consisting of many quarks and gluons, 
one will naturally ask whether it is possible to express hadronic states in terms 
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of Fock space expansion with a trivial vacuum. It is obvious that if we could 
"remove" from the theory the basic constituents with zero longitudinal light- 
front momentum, the vacuum of the full interacting theory would be the same 
as the free field theory, namely 

\vac)LF = |0) . (7) 

It must note that here "removing" from the theory the basic constituents with 
zero longitudinal light-front momentum does not mean to simply ignore dynam- 
ics of these constituents and their contributions to the bound states. Mathemat- 
ically, one can remove these constituents with zero longitudinal momentum by 
either using a prescription that requires the field variables to satisfy the antisym- 
metric boundary condition in the light-front longitudinal direction [Zhang and 
Harindranath (1993a)] or dealing with a cutoff theory that imposing a cutoff, 
fc+ > e, on the momentum expansion of each field variable, where e is a small 
number [Wilson ct al. (1994)]. Thus, the positivity of longitudinal momentum 
with such a prescription or an explicit cutoff ensures that the light-front vacuum 
must be trivial. Now a relativistic bound state can be expressed as an ordinary 
Fock state expansion: 

\M^)=f{a\h\d)m. (8) 

For QCD, a\b'^ and S are the gluon, quark and antiquark creation operators 
with nonzero longitudinal momentum, and f{a\h\S) must also be a color 
singlet operator as a polynomial function of { a\h\S }. 

ii). Light-Front Bound State Equation. Once the light-front vacuum becomes 
trivial and the light-front bound states for various hadrons are expanded in terms 
of the Fock space, the dynamic equation to determine these states is rather 
simple. Explicitly, a hadronic bound state labeled by a with total longitudinal 
and transverse momenta P+ and and helicity (the total spin along the 
longitudinal direction) A can be expressed as follows: 

|a,-P^,-PL,A) = / 2(2-jt)^' \n,XiP-^,XiP±+k±i,Xi}^>„/a{xi,k±i,Xi), 

n,\i 

(9) 

In Eq.(9), n represents n constituents contained in the state \n,XiP~^ ,XiP± + 
k±i,Xi), Xi is the helicity of the i-th constituent, and J denotes the integral 
over the space: 

^^Xi = l, and ^^k±i = 0, (10) 

i i 

where Xi is the fraction of the total longitudinal momentum that the i-th con- 
stituent carries, and k±i is its relative transverse momentum with respect to the 
center of mass frame: 

p+ 

Xi = , ki± = pi± - XiP± , (11) 
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with pf, pi± being the transverse and longitudinal momentum of the i-th con- 
stituent. ^n/a{xi,k±i,Xi) is the amplitude of the Fock state \n,XiP~^,XiP± + 
k±i, Xi) which satisfies the following normalization condition: 

E f ^J^\^n/aiXi,ku,X,)\' = l. (12) 

The eigenstate equation that the wave functions obey on the light-front is 
obtained from the operator Einstein equation = P+P^ — P^ = M^: 

HLF\a, P+, Pi, A) = P+, Pi, A) , (13) 

where H^p = P~ is the light-front Hamiltonian. Futhermore, since the boost 
on the light-front only depends on kinematics, boosting a bound state from one 
Lorentz frame to any other frame is quite simple, and is dynamically independent 
[Zhang (1994)]. Thus, if we found the bound state in the rest frame, we could 
completely understand the particle structure in any frame. This is not true in 
the instant form. In the instant form, the solutions in the rest frame are not 
easily boosted to other Lorentz frames due to the dynamical dependence of the 
boost transformation. Therefore, in each different Lorentz frame, one needs to 
solve the bound state equation of P° to obtain the corresponding wave functions. 
This is perhaps the reason why one has not established a reliable approach to 
construct relativistic wave functions in the instant field theory in terms of the 
Schrodinger picture. This obstacle is removed on the light-front. 

To see the explicit form of the light-front bound state equation, let us con- 
sider a meson wave function (for instance, a pion). The light- front bound state 
equation can be expressed as: 



{qq\H^nt\qq) {qq\Hint\qqg) 
{qqg\Hint\qq) 
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^QQg 
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(14) 



Of course, to exactly solve the above equation for the whole Fock space is still 
impossible. Practically, one has to truncate the Fock space to only include these 
Fock states with a small number of particles. For example, one may truncate all 
the high order Fock space sectors (approximately) from the valence constituent 
space. Then the light-front bound state equation is reduced to the light-front 
Bethe-Salpeter equation: 
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Note that in Eq.(15), 14// denotes an effective two-body interaction kernel. In 
other words, by "truncating" the Fock space to only keep the valence quark 
states, the complicated Eq.(14) is reduced to the manable Eq.(15) but the domi- 
nant contribution of higher Fock space to the bound states must be now described 
effectively by 14//. The residual effect should be manageable in the framework 
of perturbation theory. A true nonperturbative QCD solution to the hadronic 
bound states is if one were able to derive these effective interactions directly 
from QCD rather than that phenomenologically are put by hand. This will be 
discussed in the last Lecture. 

1.3 Phenomenological Hadronic Bound States on the Light-Front 

At the present time, how to solve for the bound states discussed above from 
QCD is still unclear. Hence, it may be useful to have some insights into the 
light-front behavior of the meson and baryon wave functions which have been 
constructed phenomenologically in describing hadrons. In fact, the phenomeno- 
logical light-front meson and baryon bound states have been studied extensively 
in the last few years, based on the light-front quark model or light-front wave- 
function description. The motivation of light-front quark model is to provide 
a simple relativistic constituent quark model for mesons and baryons that can 
yield a consistent description of the hadronic processes for both low and high 
momentum transfer. 

The general construction of the phenomenological wave functions is moti- 
vated by that of the non-relativistic constituent quark model. The constituent 
quark model has been very successful in the description of hadronic spectroscopy 
with a very simple structure, namely that all mesons consist of a quark and an- 
tiquark pair and the baryons are made of three constituent quarks, their wave 
functions satisfy the SU{Q) classification and Zweig's rule which suppresses par- 
ticle production in favor of rearrangement of constituents for hadrons [Close 
(1979)]. However, such a simple picture is very difficult to be understood within 
QCD, due to its nonrelativistic assumption and due to our belief that QCD vac- 
uum must be very complicated so that hadrons must contain an infinite number 
of quark-antiquark pairs and gluons. 

Light-front bomid states describe the relativistic hadronic structure with a 
nonrelativistic form. Furthermore, the simple vacuum state on the light-front 
ensures the validity of the Fock state expansion of hadronic states. With the 
assumption of existence of constituent quarks (of masses of hundreds of MeVs), 
the leading approximation to hadronic states that consist of a quark-antiquark 
pair for mesons and a three-quark cluster for baryons should be a reasonable 
starting point. More theoretical discussion for such a assumption from low energy 
QCD will be given later. 

However, it must note that there is a subtle problem in the description of 
hadronic structure in terms light-front bound states. That is, it is not easy to 
identify the light-front hadronic bound states with hadronic states which are 
commonly characterized by spin as a good quantum number. On the light-front. 
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we are unable to kinematically construct the hadronic bound states with fixed 
spin. The Ught-front bound states discussed in the last section are labeled by 
helicity rather than spin. In these calculations of the parton distribution and 

fragmentation functions, the hadronic bound states arc defined or classified in 
terms of the helicity. However, when we use the light-front bound states to 
compute the hadronic structural quantities, such as hadronic decay form factors 
and coupling constants, we must have states with a definite spin. A general 
solution to the spin problem on the light-front has not been found. However, 
phenomenologically, the helicity part of the bound states on the light-front can 
be transformed to a light-front spin part via the so-called Melosh transformation 
(which is exact only for free quark theory) such that the hadronic states may 
be projected (approximately but no necessary to be correct) from the set of 
light-front bound states labeled with helicities. Here, I list some meson and 
baryon light-front bound states that have been used to calculate various hadronic 
quantities in the past few years. 

The general form of the phenomenologic;al light-front hadronic bound states 
has a similar structure to the constituent quark model states: for meson states 
(with only the qq Fock space sector), 

Ai A2 

(16) 

and for baryon states (with the three quark Fock space sector), 

ip+, p^, 553) = ^ I n ^^g^o^rn-. A,) 

Xi i 

X \xi, fci_L, Ai; a;2, fc2_L, A2; l-xi- X2,-{ki^_ + k2±), A3)(17) 

where ip^^^^ is the amplitude of the corresponding qq or three quark sector (the 
wave function of the quark model): 

= J' S^^^ki±,Xi)${xi,ki^) , (18) 

with the flavor part of the wave function which is the same as in the constituent 
quark model, and S and <P are the spin and space parts that depend on the 
dynamics. By ignoring the dynamic dependence of the spin configuration and by 
using the Melosh transformation [Melosh (1974)], 

rui + XiMo - icr ■ [n X 'kij_) 

RM{ki±,mi) = , ,^x2 , ,2 

^y {rrii + XiMoY + kf^ 

where n = (0,0,1), cr is the Pauli spin matrix, mi the i-th constituent quark 
mass, and Mq satisfies 

^o=E^^^^' (20) 
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the light-front spin wave function can be given by 

E^^'{k±,XuX2) = J2{Xi\Rii{k±,mi)\si){X2\Rij{-k±,m2)\s2) 

Sl,S2 

(l/2si,l/2s2|^^3), (21) 

for mesons; for baryons the spin part is rather comphcated for a detailed con- 
struction, see for example Ref. [Schlumpf (1993)]. The momentum part of the 
wave function may be written as 



dk 

*m(a;i,fci±) =Ar„^-^exp(-kV2u;^) (22) 
for meson, where k = {k±,kz), kz = (^x — Mq — — J^-^; and for baryons 

M^i'k,^)=^»-(rTW^' (23) 

where J\f is a. normalization constant and a; is a parameter fixed by the data. 
Other phcnomcnological light-front wave functions have also been used. 

These phenomenological light-front wave functions have been widely used to 
calculate hadronic form factors and coupling constants [Chung et al. (1988)]; 
the results look pretty good for a very broad range of momentum transfer, and 
should provide a much better description than the nonrelativistic constituent 
quark model and other phenomenological descriptions. 

Nevertheless, all these arc just some phenomenological examinations of light- 
front hadronic wave functions. The true strong interaction description of hadronic 
structure is the solution of the bound state equation, Eq.(14) or approximately 
Eq.(15), from QCD. This is the main task of the recently development of QCD 
formulated on the light-front. In the remaining lectures, I will discuss the QCD 
formulation on the light-front and then explore its application to heavy quark 
systems, based heavily on the works which have been done with my collabo- 
rators in the last few years [Zhang (1993), Zhang and Harindranath (1993a), 
Zhang and Harindranath (1993b), Harindranath and Zhang (1993), Wilson et 
al. (1994), Cheung et al. (1995), Zhang (1996)]. 



2 Canonical Light-Front QCD 

2.1 Introduction 

Light-front QCD that I am going to discuss is the theory of QCD formulated on a 
light-front surface with the light-front gauge A+ = 0. Before start the discussion 
on LFQCD, I would like to make a few remarks: First of all, I would like to 
claim that any problem of QCD that can be solved in the instant formulation 
should be undoubtedly solved on the light-front. This is not surprise at all! 
However, the importance of LFQCD is that we hope to solve the subtle problems 
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in QCD that have not been solved in the instant form, such as color confinement 
and dynamical chiral symmetry breaking problems. To reach this goal, one may 
need to have some relatively complete knowledge on the canonical formulation of 
LFQCD and from which to find the key problem associated with these subtleties. 
Hence, in this lecture, I will introduce the canonical form of light-front QCD, 
then discuss the origin of the light-front gauge singularity and the light-front 
two-component formulation of QCD which has some very special structure for 
field theory that arc only manifested on the light-front. 
The QCD Lagrangian is defined by 



where F**^ = ^^'A'' - d^A^" - ig[Ai^, A% Ai" = J^a ^a^"- is a 3 x 3 gluon field 
color matrix and the T° are the generators of the SU{3) color group: [T'^, T^] = 
^jafecyc Tr(T'^T'') = i(5ab- The field variable i/j describes quarks with three 
colors and Nf fiavors, = — igA^^ is the symmetric covariant derivative, 
and m is an Nf x A^^ diagonal quark mass matrix. The Lagrange equations of 
motion are well-known; 



The following discussion and also that of the next lecture are mainly based 
on the work in collaboration with A. Harindranath [Zhang and Harindranath 
(1993a), Zhang and Harindranath (1993b), Harindranath and Zhang (1993)]. 

2.2 Light-front (phase space) Quantization 

To formulate the QCD on the light-front, the following light-front notations will 

be adopted: The space-time coordinate is denoted by x'' = {x^,x~,x±), where 
x'^ = x° + x^ is the light-front time-like component, x^ = x° — x'^ and xY (« = 
1, 2) are respectively the light-front longitudinal and transverse components. The 
light-front derivatives arc given by 9+ = = 2tj^, and d]^ = The 

product of two four- vectors is written as a • 6 = ^(a+b^ + a~b'^) — a± ■ bj_. 

The canonical theory of QCD on the light-front is constructed with the choice 
of the light- front gauge ^4+ = 0. The first question you may ask is why we choose 
the light-front gauge. The answer is as follows: 

On the light-front, the quark (more generally the fermion) fields can be de- 
coupled into i^{x) — ^p+{x) + tjj-{x) with ij)±{x) = \'y^'y'^tl){x). Then the Dirac 
equation (25) can be separated into: 



(24) 



[i-f^d^ -m + g-f^,A^)^ = , 

d^FH" + gr^'^A.^Fi^" + s^T'^r^ = , 



(25) 
(26) 



{id- + gT"A-'^il)+ = {ia± ■ £>_l + prnj^j^- , 
{id+ + gT°-A-^i}_ = {iai_ ■ £>_l + I3m}ji}+ , 



(27) 



(28) 
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where a± = 7°7_l, /? = 7°. It shows that the component V- is a constraint field 
variable, which can be solved nonperturbatively from the about equation ONLY 
if we take A+ = 0! 



V'- 



= {ia± ■ D± + Pmj t/>+ . (29) 



Secondly, due to gauge symmetry among the four components of the vector 
gauge filed, only two of them are the physically independent variables. By taking 
j4+ = 0, the equation of motions (26) for i/ = + is simply reduced to 

^(d+yA- =d+d'Ai+gpa, (30) 

where pa = f°'^'^A\d'^ A], + 2'0^T"'0_|_ is the light-front color charge density. This 
is indeed the light-front Gauss Law which can be used to determines A~ in terms 
of V'-i- and AJ,: 



)(9M;)+(— ) pX (31) 



where the operator will be defined later. It shows that with the light-front 
gauge, we can explicitly eliminate all the unphysical gauge degrees of freedom. 

Now, we can write a simple close form for the LFQCD Lagrangian in terms 
of the pure physical degrees of freedom, and (i = 1, 2), with the choice of 
the light-front gauge: 

Cqcd = \{d+Ai){d-Ai) + i^l:ld-^+ - H, (32) 

where H is the LFQCD Hamiltonian density: 

n = jid'Aif + gr^^AlAid^Ai + ^r'-^r^'^'AlAiAl^Ai 

^l{a± ■ {id± + gA_L) + /3m} W-L ' (^^J- + 9A±) + Pm}^+ 

Next we discuss the light-front quantization. A self-consistent canonical quan- 
tization requires that the resulting Hamiltonian must generate the correct equa- 
tions of motion for the physical degrees of freedom (yl^, V'+)- To reproduce 
the Lagrangian equations of motion, we need to find consistent commutators 
for physical field variables. In the light-front gauge, the LFQCD phase space 
is spanned by the field variables, AJ, , i/)-). , 1/'+ and their canonical momenta, 
£1 = l^+AJ,, 7r^_^ = fV'+i %t = —^ip+. The phase space structure which 



+ 
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determines the Poisson brackets of its variables can be found by the Lagrangian 
one- form Cdx'^ (apart from a total light-front time derivative), 

Cdx+ = ^2{SidAi + 7rv,+ #+ + d^^n^i^ 

—A\d£1 — dTr^_^_tl)+ — '^+d7r^t ) — Hdx'^ 
= ^cfr^ffdqf' - ndx+ , (34) 

where the first term on the right-hand side is called the canonical one-form of the 
phase space (note that quark fields are anticommuting c-numbers (Grassmann 
variables)). Correspondingly, the canonical equal-a;+ commutation relations are 
then given by: 

[q''{x),q^{y)U^y+=ir~^\ (35) 

Explicitly, we have 

{V+(x), ^1(2/)},+=^+ = ^A+sH^ ~ y) , (36) 

[Ai{x),Ai{y)]^^^^^=iSatS'^ (j^y^x-y). (37) 

From these commutation relations it is straightforward to verify that the Hamil- 
tonian equations of motion are consistent with Eqs.(25) and (26). As we see 
in the above light-front quantization of QCD one does not need to introduce 
the ghost field. However, this canonical formulation does not completely define 
theory for practical computations due to existence of gauge singularity. 



2.3 Light-Front Gauge Singulcirity 

The gauge singularity is perhaps the most difficult problem in non-abelian gauge 
theory that has not been completely solved since it was developed. In LFQCD, 
it arises from the elimination of the unphysical gauge degrees of freedom. To 
eliminate the unphysical degrees of freedom on the light-front, we need to solve 
the constraint equations which depend on the definition of the operator 1/9+ . 
In our formulation, we define this operator by 



e{x - Xj^)f{xj^ ,x'^,x±) , 



(38) 



where s{x) = -1, 0, 1 for .t < 0, = 0, > 0. 

In perturbation theory, the gauge singularity manifests itself clearly in mo- 
mentum space. The momentum representation of Eq.(38) is 



/(2=)=(t) I-oo^^i ■■■^^n^i^ -Xi)---e{x^_-^-xJf{xJ 



k++i€ 



+ 



(39) 
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As we see the fc+ = modes are removed with this definition. In other words, the 
singularity of is regularized. However, with such an infrared regularization, 
many infrared divergences from the small longitudinal momentum, surrounding 
the fc"*" = region, will occurs in the perturbative calculation. We will discuss 
these divergences in the next lecture. 

On the other hand, with the definition of Eq.(38), we have 

[Ai{x),Ai{y)]^^^^^ = -i^S^tS'^ e{x- - y-)S'{xi. - y^) . (40) 

This leads to the fact that A* satisfies an antisymmetric boundary condition: 

Al{x- = -oo) = -Alix- = +00). (41) 

It also shows that the zero-mode (the longitudinal momentum is zero) in is 
removed. Meanwhile, quarks in QCD should always be massive, namely their lon- 
gitudinal momentum is not really zero. Thus, with the definition of Eq.(38), the 
theory of light-front QCD does not contain zero-modes. Therefore the LFQCD 
vacuum in this formulation is always trivial! Now you may ask where is the 
nontrivial properties of QCD with such a trivial vacuum in your formulation? 

Apparently, after solving A~ component from the light-front Gauss law in 
the A+ = gauge, the gauge freedom should be completely fixed. However, a 
careful check shows that there is still a residual gauge transformation in the 
above formulation. It is given by 




C/^ = exp I - - y d'x^e''{x^)Ra{x±) |, (42) 

which is associated with the field A~ at the longitudinal infinity: 

Ra = 2^+^-1,-=^ = -j dx- [2d+d'Al{x)+gpa{x)\. (43) 

This gauge freedom can be further fixed for physical states. This is because 
the operator Ra = E~\^-^^ which is the longitudinal component of color elec- 
tric field strength at longitudinal infinity. For physical states, finite energy den- 
sity requires that the color electric field strength must vanish at the longitudi- 
nal boundary: E~\x-^j-^ = 0. This condition canonically removes the residual 
gauge freedom and leads to a constraint on the A* at longitudinal infinity: 

/oo 
dx~pa{x~,x). (44) 



The nontrivial properties of QCD in our formulation are indeed hidden in this 
condition. The main effect of this equation should be only manifested in nonper- 
turbative dynamics, i.e., in physical bound states. An explicit nontrivial eflfect 
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can be seen form the axial anomaly of QCD, for example. Consider the axial 
current equation (for zero quark mass) 

d,j^ = Nf^Tr{F^^F^n, (45) 

where the axial current is = ijj^'^^fsii'. and the dual field strength is F^"^ = 
^e^'"^''F„p. The winding number in LFQCD is defined as the net charge between 
a;+ = — oo and a;+ = oo, 

= Nf-£^J^ d^x Tr {F^^Fn . (46) 

The integration on the r.h.s. of the above equation is defined in Minkowski space 
(M) and can be replaced by a surface integral. It has been found [Zhang (1993)] 
that 

/i2 /■ x~ —CO 

AQ5 = -Nf^ dx+(fx^TT{A-[A\A^]) , (47) 

namely, a non-vanishing AQ^ is generated from the asymptotic fields of and 
their antisymmetric boundary conditions at longitudinal infinity. 

From the above canonical analysis, we can see that nontrivial features in 
LFQCD are induced by the gauge singularity and are manifested at the longi- 
tudinal infinity on the light-front. They are also associated with the light-front 
longitudinal infrared divergence in momentum space when the zero-modes are 
removed in our canonical quantization. This analysis gives us some hint where 
we should look for the problems in the study of nonperturbative QCD with a 
trivial vacuum on the light-front. 



2.4 Two-Component Formulation 

When QCD is formulated on the light-front, the theory can be expressed in terms 
of a pure two-component form. This is another useful feature of LFQCD. After 
the elimination of the unphysical gauge degrees of freedom, the QCD gauge field 
has already been reduced to the two transverse components, A\ and A^. While, 
as we will see soon that the quark field can also be written in terms of a two- 
component field (rather than the four-component field in instant form) [Zhang 
and Harindranath (1993b)]. 

To do so, we should introduce the following 7 matrix representation: 



-i 

1 



i 

1 



-le-'a-i 







(48) 



Then, one can find that the light-front project operators become 



1 
00 



00 
1 



(49) 
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V'- 







, (50) 



where (f{x) is a two-component spinor field. In the above expressions, cfi are 
the Pauli matrix. Thus, the rclativistic fermion particles can be described as 
a nonrelativistic spin 5 particle on the light-front. The canonical commutation 
relation is also reduced to: 



{fix) , ifi^y)} =5^{x-y). 



(51) 



With the above formulation, the LFQCD Hamiltonian can be rewritten as 



H 



dx (fx±{no+Tiint)=Ho + Hi, 



where 



'Hint — 'Hqqg + ^^ggg + ^qqgg + "^qqqq + ^gggg > 



(52) 

(53) 
(54) 



and 



'Hqqg = g'P^ 



+ 



1 

d+ 
1 

d+ 



1 

9+ 



qqgg 



{d-A^)+a- A^ 
a ■ d± — m)a ■ A± (p , 



{a ■ d±+ m) 



(55) 
(56) 



+2 (r^^A^a+A^) (^tya^) 



n. 



qqggi- 



K 



qqgg'2 



K 



gqqq 



1 



d+ 



1 



d+ 



(57) 
(58) 



n 



gggg 



^r".r"'I^AlAiA^Ai + 2(^) (Ald+Al) (^^) {A^d+Ai)^ 



n 



ggggi^ 



n 



gggg"^ 



(59) 



The above two-component formulation simpliiy the relativistic field theory 
structure, especially in the study of the relativistic bound state problems. 
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3 LF Time- Ordered Perturbation Theory for QCD 



3.1 About Light-Front Perturbative QCD 

Tinic-ordcrcd perturbation theory, especially the light-front timc-ordcrcd pertur- 
bation theory, provides a natural perturbative description for parton phenomena 
[Drell et al. (1970), Mueller (1989)]. The current attempts of solving nonpertur- 
bative QCD dynamics on light-front is also based on the analysis of time-ordered 
approach in Hamiltonian formulation. However, the light-front gauge singularity 
discussed above will lead to severe infrared divergences in such perturbation the- 
ory, although Eq.(39) provides a well-defined regulator (a generalized principal 
value prescription) for the small fc+ momentum. 

In covariant perturbation theory, the use of the principal value prescription 
still leads to the so-called "spurious" poles in the light-front Feynman integrals, 
which prohibit any continuation to Euclidean space (Wick rotation) and hence 
the use of standard power counting arguments for Feynman loop integrals. This 
causes difficulties in addressing rcnormalization of QCD in covariant perturba- 
tion theory with the light-front gauge. In the last decade there are many inves- 
tigations attempting to solve this problem. One excellent solution is given by 
Mandelstam and Leibbrandt, i.e., the Mandelstam-Leibbrandt (ML) prescrip- 
tion [Mandelstam (1983) and Leibbrandt (1984)], which allows continuation to 
Euclidean space and hence power counting. It has also been shown that, with 
the ML prescription, the m,ultiplicatwe rcnormalization in the two-component 
LFQCD Feynman formulation is restored [Lee and Milgram (1986)]. 

Unfortunately, the ML prescription cannot be applied to equal- a:"*" quanti- 
zation because the ML prescription is defined by a boundary condition which 
depends on a;+ itself and is not allowed in equal- a;"*" canonical theory. Yet, as we 
pointed out recently [Wilson et al. (1994)], light-front power counting differs com- 
pletely from the power counting in equal-time quantization that noncanonical 
counterterms are allowed in light-front field theory. In other words, multiplicative 
renormalization is not required in LFQCD. Furthermore, the current attempts 
to understand nonperturbative QCD in light-front coordinates is based on the 
a; "'"-ordered diagrams in which no Feynman integral is involved. Thus the power 
counting criterion for Feynman loop integrals is no longer available in LFQCD 
Hamiltonian calculations. In ^"'"-ordered perturbation theory with the principle 
value prescription, LFQCD contains severe Unear and logarithmic infrared di- 
vergences. Here I will give some results from the ^"'"-ordered perturbative loop 
calculations and renormalization of LFQCD Hamiltonian theory up to one-loop 
[Zhang and Harindranath (1993b), Harindranath and Zhang (1993)], where the 
infrared divergences are systematically analyzed. Since light-front power count- 
ing allows noncanonical counterterms, a complete understanding of renormalized 
LFQCD may not be worked out within perturbation theory; new renormalization 
and regularization approaches need to be developed, as we will see later. 
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3.2 LF x+-Ordered Perturbation Theory 

The x+-ordered perturbation theory can be obtained from the familiar pertur- 
bation expansion in quantum mechanics. The perturbation expansion of a bound 
state is given by (in the Rayleigh-Schrodinger perturbation theory): 



OO / 



^ W)] m, (60) 



Eq — Ho 



n=0 

where |<?) is a unperturbative state, Q and H'j are defined by: 

Q = |#)(>f|, H'j = Hi-AE, AE = {^\Hi\^). (61) 

With this perturbative expansion formula, the mass, the wave functions, and 
the coupling constants renormalizations can be expressed as follows. For the 
convenience of practical calculations, we consider the expressions in momentum 
space. 

i). Wavefunction renormalization: In momentum space, the perturbative 
expansion of a state is given by 



' \n,}{n,\H',\$) 



which has not been normalized, where \ni), [ri2), • • • are properly symmetrized 
(antisymmetrized) states with respect to identical bosons (fermions) in the states 
and J2 in Eq.(62) sums over all intermediate states except the initial state |^). 
The normalized wave function is defined by = ^/Z^\^), where the factor 
is the wavefunction renormalization constant: 



[P Pni ) 



-1 _ /,T,\,Tr\ _ 1 I 



ni 



ii). Mass renormalization. The mass correction can then be computed 

from the "cncrgy-lcvcl" shift, i.e., the correction to the energy of an on-mass- 
shell particle. It is obvious that the perturbative correction to the light-front 
energy (p^ ) is given by 



6p- = {^\{H - HoW) = {^Hi\^) 
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Using the mass-shell equation = p'^p~ — p5_, and recalling that p+ and p_L 
are the conserved light-front kinematical momenta, we obtain the mass renor- 
malization in the old-fashioned perturbative light-front field theory: 

= p+5p- = p+imim + p+^'^M^I^ + . . . . (65) 



iii). Coupling constant renormalization. The coupling constant renormal- 
ization is obtained by the perturbative calculation of various matrix elements 
of the vertices in Hj. Consider a vertex Hj that is proportional to the coupling 
constant g, we have 



tr Pf-Pm 
{^f\H\\m){n^\H'j\^i) 



ni 

+ E 

ni,n2 



Pi -Pni 

{$f\H'j\m){n,\H'j\n2){n,\H}]\^i) 

{Pj -Pni){pJ -Pn2 +«e) 

{^f\H'j\m){m\H}\n2){n2\H'j]\^i) 

{Pj -Pni)(pi -Pn2) 

(<?/|i/l|ni)(ni|i/||n2)(n2|i/|]|^.) 



where Zg is the multiplicative coupling constant renormalization, and Zi and Zf 
are the wavefunction renormalization constants of the initial and final states. 

It is also convenient to express the above perturbation expansion in terms of 
the diagrammatic approach. The rules for writing the expression of perturbative 
expansions from diagrams for QCD are as follows: 

• Draw all topologically distinct a;+-ordcred diagrams. 

• For each internal line, sum over hclicity and integrate using J ^^j^p^^(A;+) 
for quarks and / ^^^^^0{k+) for gluons. 

• For each vertex, include a factor of 167r^5^(p/ — Pi) and a simple matrix 
element listed in Rcf. [Zhang and Harindranath (1993b)]. 

• Include a factor {p^^ - YlmPn + [or (P/ - Y^nPn + ^^)~^] s^'Ch 
intermediate state, where Y^^Pn sum over all on- mass-shell intermediate 
particle energies. 

• Add a symmetry factor for each gluon loop coming from the sym- 
metrized boson states. 
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3.3 Perturbative Calculation of Light-Front QCD 

To illustrate the above computation scheme and to explore the severe light-front 

infrared divergences, let me list some calculations up to one-loop based on the 
j;"''-ordered diagrammatical approach. Note that besides the infrared divergence, 
which is regularized by Eq.(39), there are also ultraviolet divergences for which 
we use a transverse cut-off: yl^ > k\ > M^- Here I have also introduced a mass 
scale fj, for the minimum cut-ofF of the transverse momentum in order to avoid 
the several complicated pure infrared divergences and mass singularity from the 
massless gluon, ^ should be much larger than all other masses in the theory, and 
is considered as a renormalization scale here. 

i). Quark wavefunction and mass renormalization. The one- loop light-front 
quark energy corrections (for the three diagrams in Fig. 1, respectively) are given 

by 

-^^^ = 8^^^T+r ~' 

^P3 = -^^1 r„+i ( - 1 ) • (68) 



87r2 ^ \p+] V 2e 

This shows that, in the one- loop quark energy correction, one-gluon exchange 
gives rise to both linear and logarithmic infrared divergences. The instanta- 
neous fermion interaction contribution (see in Fig. lb) contains only one 
logarithmic divergence which cancels the logarithmic divergence in Sp'^ . The in- 
stantaneous gluon interaction contribution [Sp'^ of Fig. Ic) has a linear infrared 
divergence which precisely cancels the same divergence in 5p^ . This cancellation 
of linear infrared divergences is based on the use of the regularization for /;;+ ^ 
in Eq. (39) [Zhang and Harindranath (1993b)]. 

The quark mass correction (dropping the finite part) is then given by 

drn^ = p'^Sp~\p2^^2 = -^-^Cfm^ In , (69) 

which is longitudinally infrared divergence free; and the quark wavefunction 
renormalization constant is 



dSp 



dp 



" C,fA_21ni^V-4- (70) 



The wavefunction renormalization contains an additional type of divergence, the 
mixing of infrared and ultraviolet divergences, that does not occur in covariant 
calculations. This is the 'spurious' mixing associated with the gauge singularity. 
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Fig. 1. The x -ordered graphs for the one-loop correction to the quark mass and wave 
function renormalization. 

It corresponds to the so-called light-front double pole problem in the Feynnian 
theory with the use of the light-front gauge and the principal value prescrip- 
tion that prohibits any continuation to EucHdean space and power counting in 
Feynman loop integrals. In the a;+-ordered Hamiltonian perturbation theory the 
power counting is different. The above argument of power counting for Feyn- 
man loop integrals may be irrelevant. Furthermore, since the second order cor- 
rection to wavefunctions must be negative, the above result shows that it is the 
additional infrared divergence that gives a consistent answer for wavefunction 
renormalization. 

ii). Gluon wave Junction and mass correction. Similar calculation to the 
one-loop light-front gluon energy corrections leads to the solution: 



In the gluon sector, more severe divergences appear. It contains the quadratic 
and logarithmic UV divergences, linear and logarithmic IR divergences, and an 
unusual large longitudinal momentum logarithmic divergence. Only the linear in- 
frared divergences are cancelled with the principal value prescription. The gluon 
mass correction is not zero. The non-zero gluon mass correction of Eq.(71) is 
not surprising because it has the same divergence feature as the photon mass 
correction in light-front QED [Eq.(71) will be reduced to the photon mass cor- 
rection when we set Tf = 1, Ca = and Nf = 1]. In a. covariant calculation, 
the zero gluon mass correction is true only for dimensional regularization which 
"removes" or drops the mass correction. In the present calculation, maintaining 
zero gluon mass requires a mass counterterm, as is known in QED. The differ- 
ence between QED and QCD is only manifest in the gauge boson wavefunction 
renormalization. For wavefunction renormalization, again there is an additional 
mixing of UV and IR divergences, which again provides the correct sign for the 
wavefunction renormalization constant. 



6(1% 



1 + 




(71) 



Light-Front QCD and Heavy Quark Systems 21 

iii). Coupling constant renormalization. For convenience, we set the exter- 
nal gluon momentum q{q~^, QjJ = 0. The quark-gluon vertex is then reduced: 

Vo = 2gT^„^Sx,xA* ■ (72) 
In x^-ordered perturbation theory, the one-loop vertex correction is given by 
<5Vo = {Fi + V2 + V3 + Vi + V5 + Ve}Vo , (73) 

where Vn,n = 1 — 6 are represented the contributions from different time-ordered 
diagrams (see Ref. [Harindranath and Zhang (1993)]): 

= , ^6 = . 

To evaluate the contributions to the coupling constant we have to multiply 
Vi and V2 by ^ in order to take into account the proper correction due to 
the renormalization of initial and final states. Thus adding the contributions 
together, we have, 

S^o = (t^i + 4-^2 + V2 + V4 + V5 + Ve^ Vo 
= (^^Ca - iNfTf^ In A = ^^Vo 

Note that all mixed divergences cancel now. The correction to the coupling 
constant is given by 

9R = 9{1 + 69)=gS^l + ^ I^^Ca - jNjTj'^ In A | . (76) 

By redefining the bare coupling constant g such that ga is finite. Thus we have 
given all canonical renormalization quantities in QCD up to one-loop order based 
on the x^-ordered perturbation theory. 

Prom these results, the anomalous dimensions for quarks and gluons and the 
P function up to one-loop can be easily calculated. The anomalous dimension of 
the quark field to order is 



)ln-, (74) 



(75) 
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The momentum-dependent term implies that the quark anomalous dimension is 
gauge dependent. The anomalous dimension for the gluon field is 



which is also gauge-dependent. In the case of g"*" =0, the gauge dependent term 
can be removed, and Eq.(78) is reduced to Gross and Wilczek's result in their 
Feynman calculation with = and = [Gross and Wilczek (1974)]. The 
P function is 



which is the well-known result to one loop order and is infrared divergence free, 

as wc expected. 

From the above result, we see that there are severe light-front divergences 
in LFQCD. Systematic control of these divergences is required a priori before 
we perform any practical numerical calculation in light-front coordinates for 
QCD bound states. From the basic one-loop calculations, one can see that, in 
the a;+-ordered perturbation theory, light-front QCD involves various UV and 
IR divergences. Some of the divergences have not even been encountered in 
covariant and noncovariant Feynman calculations to the same order. Among 
various light-front divergences, there are two severe divergences one has to deal 
with in the x^-ordered theory for light-front QCD. The first is the mixing of UV 
and IR logarithmic divergences in wavefunction renormalization. The occurrence 
of the mixing divergences may not be a severe problem. The mixing divergences 
should be cancelled completely for physical quantities, as we have seen from 
the coupling constant renormalization. We expect that the problem of mixing 
divergences may not exist when we consider real physical processes. The second 
problem is the infinite gluon mass correction. In the time-ordered perturbation 
theory dimensional regularization is not available to avoid the nonzero gluon 
mass correction. To have a massless gluon in perturbation theory, we have to 
introduce a gluon mass counterterm. In the leading order (one-loop) calculation, 
there is no difficulty arising from a gluon mass counterterm. However, when we 
go to the next order, it has been found that the gluon mass counterterm leads to a 
noncancellation of infrared divergences. The non- vanishing infrared divergences 
could introduce non-local countcrtcrms in both the longitudinal and transverse 
directions. In instant quantization, such non-local counterterms are forbidden for 
a renormalizable theory. Here, these non-local counterterms are allowed by the 
light-front power counting. This is a special feature of LFQCD. One speculation 
from this property is that the non-local counterterms for infrared divergences 
may also provide a source for quark confinement [Wilson et al. (1994)]. 

In summary, renormalization in LFQCD Hamiltonian theory is very differ- 
ent from conventional Feynman theory and it is an entirely new subject where 
investigations are still in their preliminary stage. In perturbative calculations. 




(78) 




(79) 
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careful treatment could remove all severe infrared divergences for interesting 
physical quantities in LFQCD. For nonperturbative studies, the cancellation of 
severe infrared divergences may not work because certain approximations (e.g., 
Fock space truncation) might bo used. These approximations may also break 
many important symmetries such as gauge invariance and rotational invariance. 
It is the hope of the current investigation of light-front renormalization theory 
that the countertcrms for the light-front infrared divergences may restore the 
broken symmetries and also provide an effective confining LFQCD Hamiltonian 
for hadronic bound states. 

4 Light-Front Heavy Quark Effective Theory (HQET) 

4.1 About Heavy Quark Symmetry and HQET 

The rich information about electrowcak and strong interactions that can be 
extracted from various heavy hadron decays has led to the extensive exploration 
of the QCD based and modcl-indopcndont description of heavy hadrons in the 
past few years. This is mainly due to the discovery of heavy quark spin- flavor 
symmetry (HQS) in heavy meson decays by Isgur and Wise [Isgur and Wise 
(1989), Neubcrt (1994)]. For a typical example, with the HQS, all six form factors 
m B ^ D and B D* decays are reduced to an universal function, called 
the Isgur- Wise function, and the normalization of this universal function at the 
zero-recoil point provides a model-independent determination of the Kabayshi- 
Makawa matrix element \Vcb\- Similarly in heavy-baryon decays, the application 
of HQS also leads to tremendous simplifications. 

On the other hand, heavy quark symmetry can be derived from QCD in heavy 
mass limit mg oo, via the so-called heavy quark effective theory (HQET) 
[Eichten and Hill (1990), Georgi (1990)]. The later is an effective theory of QCD 
for heavy quark expansed in inverse powers of heavy quark mass mg. In fact, 
HQET provides us with a systematical expansion of QCD dynamics in terms of 
the dimensionless parameter Aqcd /niq. and it serves as a theoretical framework 
for the systematical computation of the l/mg corrections to the limit mg — > oo. 
Thus, the HQET offers us a new channel to explore the intrinsic properties of 
hadronic structure from QCD. 

In order to actually compute any physical observables and make definite 
predictions, one still has to confront the non-perturbative QCD dynamics. Cur- 
rently, except for the lattice approach, the main physical quantities, such as 
Isgur- Wise function, can only be computed in various hadronic models, such as 
the constituent quark model, the bag model, and QCD sum rules. It would be 
very interesting if one could calculate the Isgur- Wise function, or any hadronic 
form factors, directly from QCD. This requires to construct explicitly the heavy 
hadron bound states within the HQET, which is also necessary for a complete un- 
derstanding of heavy hadron dynamics. We are motivated by such requirement to 
reformulate HQET on the light-front, from which we hope to consistently study 
the heavy hadron bound state problem [Cheung et al. (1995)]. Meanwhile, as we 
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know for light quark systems, both quark confinement and spontaneously chiral 
symmetry breaking play an essential role to the quark dynamics in hadrons. In 
order to provide a nonperturbative QCD description for light quark systems, 
it is necessary to understand the underlying mechanism for quark confinement 
as well as for chiral symmetry breaking. This will certainly make the problem 
most complicated. However, for heavy quark systems, chiral symmetry is explic- 
itly broken so that confinement is the sole nontrivial feature influencing heavy 
quark dynamics. Choosing the heavy hadron systems should be a good starting 
point in the study of nonperturbative QCD. The light-front HQET discussed 
here is mainly based on the works collaborated with C. Y. Cheung and G. L. 
Lin [Cheung et al. (1995)]. 

4.2 1/mQ Expansion of the Heavy Quark Lagrangian on the LF 

Let us begin with the QCD Lagrangian for a heavy quark: 

£ = Qiijp - mQ)Q, (80) 

where Q is the heavy quark field operator, mg the heavy quark mass and 
the QCD covariant derivative. 

In the instant formalism, HQET is obtained by redefining the heavy quark 
field as: 

Q{x) = e-'^^---[h.,{x) + H,{x)l (81) 

where v is the four velocity of the heavy quark, such that v"^ = 1; hy{x) and 
Hy{x) are respectively the so-called large and small components of the heavy 
quark field, satisfying /jhj,{x) = ft„(.x) and ^Hy{x) = —Hy{x). From the QCD 
equation of motion, one can express Hy(x) in terms of hy{x) and show that 
the former is suppressed by l/mq compared to the later. Using Eq.(81) and 
the relation between hy{x) and Hy{x), one can systematically expand the QCD 
Lagrangian in powers of l/mg, and arrive at an effective theory for the heavy 
quark. 

In the framework of light-front quantization, the situation is quite different. 
Before taking the heavy quark mass limit, the quark field is already divided into 
two parts: Q{x) = Q+{x) + Q-{x), with Q±{x) = A±Q{x) = i7°7±Q(a;). The 
Dirac equation for Q can then be rewritten as two coupled equations for Q±: 

iD-Q+{x) = {ia^_ ■ + (3mQ)Q_{x), (82) 
iD+Q_ (x) = (ia_L • £>_L + pmQ)Q+ {x) , (83) 

where a± = "f^jx and /? = 7°. As we known only the plus-component Q+{x) 
is the dynamical field. The minus-component Q-{x) is a light-front constraint 
that can be determined from Q+{x). In terms of Q+{x), the QCD Lagrangian 
(1) for the heavy quark can be rewritten as 

£ = QliD-Q+ - Ql{ia± ■ D± + PmQ)Q-, (84) 
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where Q- can be eliminated by Eq.(83). 

To derive the hght-front HQET, we use the same redefinition of the heavy 
quark field as in the covariant case, 

Q(a;) =e-^'"«'""Q,(a;), (85) 

but without imposing any constraint on the new variable Q„ to separate the large 
and small components. It follows that Q±{x) = e~""'2^'^Q„±(a;). Substituting 
this result into Eq.(83), we obtain 

Qv-{x) = _^ ^ ia± ■ D±_ + mQ{a±_ • w± + /3) Qv+{x). (86) 

mQV+ + iU+ L J 

It is worth noting that in the ordinary light-front formulation of field theory, 

the elimination of the dependent component requires the choice of the light- 
front gauge = 0, and a specification of the operator 1/9+ which leads to 
severe light-front infrared problem that has still not been completely understood 
[Zhang and Harindranath (1993a)]. However, for the heavy quark field with the 
redefinition of Eq.(85), the above problem does not occur since the elimination 
of the dependent component Qy- now depends on the operator 1/ {mQV^ +iD'^) 
which has no infrared problem. Moreover, it has a well defined series expansion 
in powers of iD^ /niQ: 

1 _ 1 V f 1 VY .D+y-^ 

mQV++iD+ v+ ^,\mnJ \ v+ J ' 

n=l 

Thus, the heavy quark QCD Lagrangian (80) can be expressed in terms of Qy+ 
alone. The complete l/mg expansion is given by 



CSO 

Co + J2^n, (88) 



where 



a-D = a^-JP^- "^-^^+V . (89) 



This is the light-front efi'ective heavy quark Lagrangian. 
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4.3 Properties of the Light-Front HQET 

In the symmetry limit, the Ught-front HQET reduces to 

Co = ^Ql+{iv-D)Q,+ , (90) 

which clearly exhibits the flavor and spin symmetries, because it is independent 
of Dirac 7-matrices and the heavy quark mass, as in the covariant formulation. 

However, beyond the symmetry limit, the light-front HQET has several ad- 
vantages over the instant formulation. In the instant HQET, the non-leading 
terms contain high order time-derivatives; consequently it is difficult to per- 
form a consistent canonical quantization beyond the limit mq 00 [Lcbcd and 
Suzuki (1991)]. It is remarkable to see that in the light-front HQET, only linear 
time-derivative appears, and it resides in Cq. The presence of the matrix fi in 
the non-leading terms eliminates all light-front time derivative terms. This can 
be seen more clearly in Eq.(88). Thus the canonical quantization of light-front 
HQET is straightforward: First of all, the canonical conjugate of the dynamical 
variable Q„+ is given by 

dC 

which docs not involve any l/mq corrections. Then using the light- front phase 
space quantization [Zhang and Harindranath (1993a)], we obtain the basic anti- 
commutation relation: 

{Qv+{x) , Ql+{y)},+^y^ = A+6Hx - y), (92) 

which is valid to all orders in I/ttiq. 

The second very useful property of the light-front HQET is that the heavy 
quark effective Hamiltonian is well defined on the light-front. From Eqs.(88) and 
(91), we obtain the light-front heavy quark effective Hamiltonian, 

H = j dx-(fx±n{x) (93) 

with the Hamiltonian density 7i given by 

H = ■:^Ql+{v-d+ - 2v^ ■ 5x)Q.+ - ^Ql+{v ■ A)Q,+ +Hm^ (94) 

and 

00 00 

71=1 n—1 

This light-front heavy quark effective Hamiltonian can serve as a basis for con- 
structing heavy hadron bound states, as we will see in the next lecture. It is 
also useful for the study of the l/mg corrections in heavy quark dynamics. 
Specifically, suppose we choose the light-front gauge {A'^ = 0) in the light-front 
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HQET, we see immediately from Eq.(88) that, in the symmetry breaking terms, 
the power of the gluon field does not increase with that of l/mg. This prop- 
erty, which is unique to the light-front formulation, may greatly simplify our 
treatment of 1/ttiq corrections. Meanwhile, note that the non-leading light- front 
effective Hamiltonian 7i„ is precisely the minus of the corresponding effective 
Lagrangian £„ given by Eq.(88). This simple relation is not valid in the instant 
HQET, due to appearance of the high-order time-derivative terms. 

Furthermore, since we have not chosen any specific gauge, and also there is no 
light-front infrared divergent problem for the heavy quark sector, short-distance 
QCD corrections to the heavy quark current and the effective Lagrangian must 
be the same as those calculated in the covariant formulation. Of course, an 
explicit calculation of the short-distance effects in the light-front HQET is needed 
to confirm the above statement, which has not been done as I known. 

4.4 Isgur-Wise Function 

The heavy quark current can also be systematically expanded in l/mq on the 
light-front. In the heavy mass limit, it reduces to the following familiar from: 

Q^{x)rQ'{x) = e-'^^Q^^'-^Q^^^^/^f (a;)r/if (a;), (96) 

where = |l + "^'y+'^^ ^Qv+- Consequences of the spin symmetry can be 
readily derived using this zeroth order heavy quark current. As an example, 
consider the matrix elements 

{PQ,iv')\!Crhf\PQ.{v)) and {P^,{v')\hi';rK''\PQ4v)), (97) 

where F stands for any arbitrary gamma matrix, Pq and Pq represent respec- 
tively a pseudoscalar meson and a vector meson containing a single heavy quark 
Q. The quantum numbers of the heavy mesons can be efficiently accounted for 

by the interpolating fields: \PQi{v)) = h], 75^^10), \PQi{v)) = h\ ^ ^v\Q), where 
e is the polarization vector of the vector meson, and represents the fully in- 
teracting hght quark (or brown muck). Prom (0|Q„+Q|^_|_|0) = ■^^+, it is easy 
to show that 

.fi;.(l + 2ij;i±i!)^.,(l + £ij;i±i!)„.i±^^ (98) 

Hence, in the heavy mass limit, the heavy meson decay matrix elements on the 
light-front take the familiar forms: 

{PQ,{v')\hi':rK'^\PQ.{v)) 
{p*Q,W)\t'rhf\PQ.{v)) 



= Tr{,.{^)r[^),.M} (99) 
= T.{^(l^)r(i±^),3M}. (100) 
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where M is the transition matrix element for the Hght quark [Wise (1991)], 

M = (0|4'4|0) ^{v' ■ v)I. (101) 

Thus spin symmetry impHes that the transition matrix elements (97) are de- 
scribed by a single form factor £^{v ■ v'), which is just the famous Isgur-Wise func- 
tion. An explicit calculation of the Isgur-Wise function from light-front bound 
state wave function will be given in the next lecture. 

5 Quark Confinement and Heavy Hadron Bound States 

5.1 A Weak-Coupling Treatment to Nonperturbative QCD 

There are two fundamental problems in QCD for hadronic physics, the quark 
confinement and the spontaneous breaking of chiral symmetry. These two prob- 
lems are the basis for solving the low-energy hadronic bound states from QCD 
but none of them has been completely understood. Recently, Wilson et al. pro- 
posed a now approach to determine hadronic bound states from nonperturbative 
QCD on the light-front with a weak-coupling treatment (WCT) [Wilson et al. 
(1994)]. The key to eliminating necessarily nonperturbative effects is to con- 
struct an effective QCD Hamiltonian in which quarks and gluons have nonzero 
constituent masses rather than the zero masses of the current picture. The use 
of constituent masses cuts off the growth of the running coupling constant and 
makes it conceivable that the running coupling never leaves the perturbative 
domain. The WCT approach potentially reconciles the simplicity of the con- 
stituent quark model with the complexities of QCD. The penalty for achieving 
this weak-coupling picture is the necessity of formulating the problem in light- 
front coordinates and of dealing with the complexities of renormalization. 

Succinctly, this new approach of achieving a QCD description of hadronic 
bound states can be summarized as follows: Using a new renormalization scheme, 
called similarity renormalization group (SRC) scheme that is recently proposed 
by Glazek and Wilson [Glazek and Wilson (1994), Wilson et al. (1994)], one can 
obtain an effective QCD Hamiltonian H\ which is a series of expansion in terms 
of the QCD coupling constant, where A is a low energy scale. Then one may solve 
from Hx the strongly interacting bound states as a weak-coupling problem. The 
WCT scheme contains the following steps: (i) Compute explicitly from SRG the 
Hx up to the second order and denote it by Hxo as a nonperturbative part 
of Hx- The remaining higher order contributions in Hx are considered as a 
perturbative part Hxi- (ii) Introduce a constituent picture which allows one to 
start the hadronic bound states with the valence constituent Fock space. The 
constituent quarks and gluons have masses of a few hundreds MeV, and these 
masses are functions of the scale A that must vanish when the effective theory 
goes back to the high energy region, (iii) Solve hadronic bound states with Hxo 
nonperturbatively in the constituent picture and determine the scale dependence 
of the constituent masses and the coupling constant. The coupling constant g 
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now becomes an effective one, denoted by g\. If we could show that with a 
suitable choice of A at the hadronic mass scale, the effective coupling constant 
g\ can be arbitrarily small, then WCT could be applied to H\ such that the 
corrections from Hxj can be truly computed pcrturbativcly. If everything listed 
above works well, we may arrive at a weak-coupling QCD theory of the strong 
interaction for hadronic bound states. 

With the idea of SRG and the concept of coupling coherence; [Pc;rry and 
Wilson (1993)], Perry has shown that upon a calculation to the second order, 
there exists a logarithmic confining potential in the resulting LFQCD effective 
Hamiltonian [Perry (1994)]. This is a crucial finding to light-front nonperturba- 
tive QCD. However, the general strategy of solving hadrons through the WCT 
scheme is far to be completed. Very recently, I used SRG to analytically de- 
rive from the light-front HQET [Cheung ct al. (1995)] a heavy quark QCD 
Hamiltonian which is responsible to heavy hadron bound states. The resulting 
Hamiltonian explicitly contains a confining interaction between a heavy quark 
and a heavy antiquark at long distance plus a Coulomb-type interaction at short 
distance. With this effective QCD Hamiltonian, I study the strongly interact- 
ing heavy hadronic bound states, from which I can provide a WCT to nonper- 
turbative QCD on the light- front, at least for heavy quarkonia. The following 
discussion is mainly based on my recent work, Ref. [Zhang (1996)]. 



5.2 Light-front SimilEirity Renormalization Scheme 

The basic idea of the SRG approach is to develop a sequence of infinitesimal 
unitary transformations that transform an initial bare Hamiltonian to an 
effective Hamiltonian H\ in a band-diagonal form relative to an arbitrarily cho- 
sen energy scale A: 

Hx = SxH^'sl (102) 

Here the band-diagonal form means that the matrix elements of H\ involving 

energy jumps much larger than A will all be zero, while matrix elements involving 
smaller jumps or two nearby energies remain in Hx. The similarity transforma- 
tion should satisfy the condition that for A — > oo, H\ — > and S\ ^ 1. 

Here, I shall follow the formulation of SRG developed on the light-front [Wil- 
son et al. (1994)]. The effective Hamiltonian we seek is Hx with A being of order 
a hadronic mass (^ 1 GeV). We begin with a given bare Hamiltonian which can 
be written by = Hq + Hf, where Hq is a bare free Hamiltonian and Ei is its 
eigenvalue. Consider an infinitesimal transformation, then Eq.(102) is reduced 
to 

^ = [Hx,T,], (103) 

which is subject to the boundary condition limx^oo H\ = . 

To force the Hamiltonian H\ becoming a band-diagonal form in energy space, 
we need to specify the action of the generator operator Tx. This can be done by 
introducing the scale A with xxij — ^Ijf^lx ^"^^^ ^ smearing function fxij = 
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f{x\ij) such that when a; < 1/3, / = 1; when x > 2/3, / = 0; and / may be a 
smooth function from 1 to for l/3<x<2/3. We can write H\ = Hq + Hj\ 
because Hq is invariant under transformations. Then Eq.(103) can be reexpressed 
as 

= fxij[Hix,Tx]ij + —i\nfxij)Hxij, 

= E~e; {^^~ ^*^) [^^A, Txh - ^ (In fxij)Hxij | . (104) 

Here we have used the notation Aij — {i\A\j). and |?') is an cigcnstatc oi Hq. Since 
f{x) vanishes when x > 2/3, one can see tliat Hxij does indeed vanish in the 
far off-diagonal region. It also can be seen that Txij is zero in the near-diagonal 
region. The solutions for Hjx and T\ are 

Hix = Hf^ + [Hix',Tx'] , Tx = Hf^T + [HiX',Tx'] , (105) 
^ V 'r ^ . 't 

where Hf,,. = h^.Hf,., Hf^^,. = (A/,..) Hf,., and 

/•OO 

Xx'ii =-fxij / d\'Xy^,, (106) 

Finally, one obtains an iterated solution for Hx, 
Hx = (ho + Hf)j + [[Hfy^Hfyj]^ 

+ [[Hfx" ) Hfy,rp\ , Hfyj] + [Hfy , [Hfy, , Hfy,j] ] + . . . 

V ' « 'R' ' « 'T' I 

^ « 'r « 'r 

= iff ^ + +hP + ..., (108) 

Thus, through SRG, we eliminate the interactions between the states well- 
separated in energy and generate the effective Hamiltonian of eq.(108). The 
expansion of eq.(108) in terms of the interaction coupling constant brings in or- 
der by order the full theory corrections to this band diagonal low energy Hamil- 
tonian. 

Explicitly, the bare Hamiltonian input in the above formulation can be 
obtained from the canonical Lagrangian with a high energy cutoff that removes 
the usual UV divergences. For LFQCD dynamics, the bare Hamiltonian has been 
constructed in lecture H (for detailed discussion, see [Zhang and Harindranath 
(1993a), Zhang and Harindranath (1993b)]). Instead of the cutoff on the field 
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operators which is introduced in ref. [Wilson et al. (1994)], I shall use a vertex 
cutoff to every vertex in the bare Hamiltonian [Zhang (1996)]: 

e{AyP+-\pr-pj\), (109) 

where p~ and pj are the initial and final state light-front energies respectively 
between the vertex, A is the UV cutoff parameter, and P+ the total light-front 
longitudinal momentum of the system we arc interested in. Eq.(109) is also 
called the local cutoff in light-front perturbativc QCD [Lepage and Brodsky 
(1980)]. All the yl-dependences in the final bare Hamiltonian are removed by the 
counterterms. The use of eq.(109) largely simplifies the analysis on the cutoff 
scheme in ref. [Wilson et al. (1994)]. 

Meanwhile, in SRG calculation, we should also give an explicit form of the 
smearing function f\ij. One of the simplest smearing functions that satisfies the 
requirements of SRG is a theta- function [Zhang (1996)]: 

fxij=0{l-xxij). (110) 



On the light- front, it is convenient to redefine x\ji = — '— — — -. Then we 
can further replace the above smearing function by the following form: 

fxij = o{§^-\^p^M (111) 

where '^P~j = — P^ is the light-front free energy difference between the 
initial and final states of the physical processes. The light-front free energies of 
the initial and final states are defined as sums over the light-front free energies 
of the constituents in the states. 

With the definition of (111), Eq.(108) can be reduced to 

Hxii = - \AP-.\) [nf, +Y.Hp,,Hg, + ^1 + . . .1 (112) 



The front factor (the theta-function) in the above equation indicates that Hx 
only describes long distance interactions (with respect to the scale A) which is 
responsible to hadronic bound states. The function gxij in eq.(112) is given by 

gXijk = / d{X' /P^)f\'ik ,/,,2 / p. ^ A'jfe 

JX'/p+ d{\' /P+) 

= e{\APj,\ - \^/P+)6{\AP7^\ - \APr^\). (113) 
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5.3 Heavy Quark Confining Interaction 

Now we can use SRG to the light-front HQET to derive a heavy quark confining 
Hamiltonian, from which we may solve from QCD the heavy hadron bound states 
directly. 

In the large mq limit, only the leading (spin and mass independent) Hamilto- 
nian is remained. The l/rrig terms (n > 1) in (88) can be regarded as perturba- 
tive corrections to the leading order operators and states. To determine confining 
interactions in heavy quark systems, the leading heavy quark Hamiltonian plays 
an essential role. With the light-front gauge = 0, the leading-order bare QCD 
Hamiltonian density is 



where V'+ is either the heavy antiquark field or the light-front quark field operator 
in the present consideration. Note that besides the leading term in eq.(88), the 
above bare Hamiltonian has also already included the relevant terms from the 
gauge field part, — |Tr(F^^F'^''), of the QCD Lagrangian. These terms come 
from the elimination of the unphysical gauge degrees of freedom, the longitudinal 
component A~ [Zhang and Harindranath (1993b)]. Eq.(114) has obviously the 
spin and flavour heavy quark symmetry, or simply the heavy quark symmetry. 

The above leading Hamiltonian (or Lagrangian) is the basis of the QCD- 
based description for heavy hadrons containing a single heavy quark, such as B 
and D mesons. As recently pointed out by Manuel et al. [Manuel and Schuler 
(1995)] the purely heavy quark leading Lagrangian may be not appropriate to 
describe heavy quarkonia. This is because the anomalous dimension of QCD 
radiative correction to QQ currents contains an infrared singularity in the limit 
of two heavy constituents having equal velocity. Such an infrared singularity is a 
long distance effect and should be absorbed into quarkonium states. To avoid this 
problem, they argued that one may incorporate the effective Hamiltonian with 
at least the first order kinetic energy term into the leading Hamiltonian [Manuel 
and Schuler (1995)]. The light-front kinetic energy can be obtained from eq.(95). 



As a consequence, in the heavy mass limit, quarkonia have spin symmetry but 
no flavour symmetry. 

i). Confining Hamiltonian for Heavy Quarkonia. Within light-front HQET, 
we now follow the procedure described above to find an effective QCD Hamil- 
tonian for QQ systems. The bare Hamiltonian for QQ systems contains (114) 



Hid = ^Ql+{v-d+ - 2v^ ■ dA_)Q. 
1,11+ ^ 




(114) 




(115) 
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and (115) for both heavy quark and antiquark plus the full QCD Hamiltonian 
for gluons and light quarks [Zhang and Harindranath (1993b)]. Since the kinetic 
energy (115) should be the same order of the Coulomb interaction, we may treat 
the kinetic energy in the same way as the instantaneous QQ interaction [the last 
term in eq.(114)]. Thus, the free Hamiltonian Hq used in SRG is given only by 
the first term in eq.(114) plus the free gluon Hamiltonian. 

With the above consideration, it is easy to compute the effective Hamiltonian 
eq.(112) for QQ systems. Following the WCT ideas, we shall calculate Hx for 
QQ systems up to the second order in the initial and final states defined by 
\i) = 6j(fci, Ai)rft (fc2, A2)|0) and \j) = 4(^3, A3)dJ(A:4, A4)|0), respectively, where 
ki is the residual momentum of heavy quarks, = rnqv^ + k^, and Ai its 
helicity. The result is 

where 



A 



2ni + A\2y' - 2t/ + 1) -i - 2^Cf^lne I, (117) 



^QQij ~ '^-L - '^'-l) = 2(27r)^53(fci + ^2 - fcs - ki)Sx,xJx^x. 

^ -V(r-)(r'^) 



a' 



+ =2^{A{y -y',K^-K'^,A)-X^)}. (118) 

(Ki-02 + (y-t/')2/ J 

Here we have introduced the longitudinal residual momentum fractions and the 
relative transverse residual momenta, 

y = k^/K+ , K± = ki± - yK±, 

y' = k+/K+ , k1 = k3± - y'K^, (119) 

where K^^ is defined as the residual center mass momentum of the heavy quarko- 
nia: = Av^, and A = Mh — tuq — rn-g is a residual heavy hadron mass. It 
follows that = fc+ + fc+ = fc.+ + , = ki± + k2± = fcsx + k4±. Since 
< p'^ = m.Qi'^ + < MffV^, in the heavy quark mass limit, we have 
Mh 2m,Q so that —rnqv^ < k^ , k^ < rnqv^ . Hence, the range of y and y' 
are given by — oo < y,y' < oo. We have also defined in eq.(118) 

A{y - y', Ac^ - n'^,A) ^ ^'^^~_''y^^ + 1^ " y'\^^- (^^O) 

Eq.(116) is the nonperturbative part of the effective Hamiltonian for heavy 
quarkonia in the WCT scheme, in which we have already let UV cutoff param- 
eter yl ^ oo and the associated divergence has been put in the mass correction. 
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The kinetic energy (115) is now included in the above effective Hamiltonian 
[the l/mg term in eq.(117)]. Note that there is an infrared divergent term in 
eq.(117) which comes from the quark self-energy correction in SRG, where e is 
an infrared cutoff of the momentum fraction / ^ and the; longitudinal 
momentum carried by gluon in the quark self-energy loop. The usual mass cor- 
rection 5mQ — -^CfA In has been rcnormalizcd away in eq.(117). In the 
WCT scheme, by removing away this mass correction, we should assign the cor- 
responding constituent quark mass in iJ^o being A-dependent. But, the heavy 
quark mass is larger than the low energy scale. Its dependence on A should be 
very weak and could be neglected. While, the QQ interaction (118) contains 
two contributions: the instantaneous interaction plus the second order contribu- 
tion in eq.(112) [i.e. the terms proportional to the theta function in eq.(118)]. 
We shall show next that the above Vq-q is indeed a combination of a confining 
interaction plus a Coulomb-typc interaction. 

ii). Quark Confinement on the Light-Front. In our framework, LFQCD vac- 
uum is trivial. The nature of nontrivial QCD vacuum structure, the confinement 
as well as the chiral symmetry breaking, must made manifestly in H\ in terms 
of new effective interactions. We will see that Hxo explicitly contains a confining 
interaction at long distances. The interactions associated with the chiral symme- 
try breaking may be manifested in the fourth order computation of H\ for light 
quark systems [Wilson et al. (1994)], but these interactions are not important 
in the study of heavy hadrons here. 

The confining interaction can be easily obtained by applying the Fourier 
transformation to the first term in (118). It is convenient to perform the calcu- 
lation in the frame K± = 0, in which 



f 

Jo 



27r2 



q±m 2Ji(|x_L|g±TO) 



\Xj_\q±r, 



(121) 



where we have used the relation q'^ = — = K^{y — y'), q± = ki± — k3± = 

Kj_ — k'j_ for K± = 0, while qj_jn = \J 'W^l^ ~ TT^I^^ ^^'^ Ji{x) is a Bessel 
function. An analytic solution to the integral (121) may be difficult to carry out. 
However, the nature of confining interactions is a large distance QCD behavior. 
We may consider the integral for large and x^. In this case, if q'^x~ and/or 

|a;_L|9_Lm are large, the integration vanishes, yet J\{x) = f + fg H for small x. 

The dominant contribution of the integral (121) for large x~ and x±_ comes from 
the small q'^ such that g+x" and/or |x_l|(J'_lto must remain small, which leads 
gi9+x- ■2J^(\x^\q^_m) ^ ^ Ti^ig corresponds to g+ < and/or q+ < j-Kj^. 
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K+ 



\{q+ <^ < i^Tpp-, eq.(121) is reduced to 



9x 
27r2 



where a term ~ is neglected since x~ is large, and e is an infrared cutoff of 
the momentum fraction jK^ . It is the same as the divergence occurs in the 
quark self-energy contribution so that the above infrared logarithmic divergence 
(~ Ine) exactly cancels the divergence in eq.(117) for color single states. What 
remains is a logarithmic confining interaction (except for a color factor): 

Vooniix- ,x^) ^ ^^|^ln(X+|x-|). (123) 

Similarly, when g+ < ^Jf^^^^i < we have 

(124) 

where the term ~ has also been ignored because of the large a;^ . Again, the 

infrared divergence (~ Ine) is cancelled in Hx for physical states, and we obtain 
the following confining interaction: 

V,onfXx-,x^) ^ ^^ln(A2|xxn. (125) 

Hence, the effective Hamiltonian Hxq exhibits a logarithmic confining interaction 
between a heavy quark and a heavy antiquark in all the directions of a;~ and x± 
space. 

The Coulomb interaction corresponds the second term in (118), its Fourier 
transformation (except for the color factor) is 

^ ^ J_ /■ da;-d2a:^e^(--«"+«---) ( ^] (126) 



where ri = \jx\+ {x~Y which is defined as a "radial" variable in the 

light-front space [Wilson ct al. (1994)]. Eq.(126) shows that the Coulomb inter- 
action on the light-front has the form 

Fco„r(.-,.^)~-|iAi. (127) 

Thus, we have explicitly shown that Hxq contains a Coulomb interaction at short 
distances and a confining interaction at long distances. 

Moreover, a clear light-front picture of quark confinement emerges here. To be 
specific, we define quark confinement as follows: i) There is a confining interaction 
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between quarks such that quarks cannot be weh-separated; ii) No color non- 
singlet bound states exist in nature, only color singlet states with finite masses 
can be produced and observed; and iii) The conclusions of i-ii) are only true for 
QCD but not for QED. 

We have shown explicitly the existence of a confining interaction in Hxq. One 
can also easily see from H\q the non-existence of color non-singlet bound states. 
This is essentially related to the infrared divergences in H\q. From eqs.(122) and 
(124), we find that the uncancelled instantaneous interaction contains a logarith- 
mic infrared divergence. Except for the color factor, this infrared divergence has 
the same form as the divergence in eq.(117). Thus, we immediately obtain the 
following conclusions. 

(a) . For a single (constituent) quark state, the interaction part of H^q does 
not contribute to its energy. The remaining infrared divergence from quark self- 
energy correction implies that the dynamical quark mass for a single quark state 
is infinite (infrared divergent) and cannot be renormalized away in the spirit of 
gauge invariance. Equivalcntly speaking, single quark states carry an infinitely 
large mass and therefore they cannot be produced. 

(b) . For color non-singlet composite states, the color factor (T")q^(T")5^ in 
the QQ interaction is different from the color factor C/ = Tr(T"T"). Therefore, 
the infrared divergence in the self-energy correction also cannot be cancelled by 
the corresponding divergence from the uncancelled instantaneous interaction. As 
a result, color non-singlet composite states are infinitely heavy that they cannot 
be produced as well. 

(c) . For color singlet QQ states, the color factor (T")(r'') Cf. Thus, the 
infrared divergences are completely cancelled and the resulting effective Hamil- 
tonian is finite. In other words, only color singlet composite are physically ob- 
servable. 

Finally, we argue that the above mechanism of quark c;onfinement is indeed 
only true for QCD. As we have seen the light-front confinement interaction is 
just an effect of the non-cancellation between instantaneous interaction and one 
transverse gluon interaction generated in SRG. Such a non-cancellation arises in 
SRG because we introduce the energy scale A. Introducing the energy scale A in 
SRG forces the transverse gluon energy involved in the QQ effective interaction 
never be less than a certain value (the energy scale A). This implies that the 
gluon may become massive at the hadronic mass scale. Of course, such a gluon 
mass must be a dynamical mass generated from the highly nonlinear gluon in- 
teractions. In other words, the above confining picture is indeed a dynamical 
consequence of non-Abelian gauge theory. This confinement mechanism is not 
valid in QED. In QED, since photon mass is always zero, the photon energy cov- 
ers the entire range from zero to infinity. Thus, in QED, we can always choose 
the energy scale A being zero. With A = 0, the infrared divergences do not occur 
in the electron self-energy correction. As a result, the renormalized single elec- 
tron mass is finite, in contrast to the divergent mass of single quark states. For 
the same reason, with A = 0, the instantaneous interaction in the effective QED 
Hamiltonian is also exactly cancelled by the same interaction from one trans- 
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verse photon exchange so that only one photon exchange Coulomb interaction 
remains. Thus, applying SRG to QED and let A = in the end of procedure, 
we obtain a conventional QED Hamiltonian which only contains the Coulomb 
interaction. 

iii). Extension to Heavy- Light Quark Systems. We can also apply SRG to 
the heavy-light quark system (heavy hadrons containing one heavy quark). The 
bare cutoff Hamiltonian we begin with for heavy- light quark systems is the com- 
bination of the heavy quark effective Hamiltonian (114) and the full Hamilto- 
nian for the light quarks and gluons. We may also introduce the residual center 
mass momentum for heavy- light systems, = Av~^ = pf + kf = pt f^t ^ 
K± = Av± = pi± + ki± = p2± + k2±, where A = Mh — mq, pi and p2 are 
the light antiquark momenta and ki and fc2 the residual momenta of the heavy 
quarks in the initial and final Qq states respectively. 

Following the general procedure, it is easy to find the nonperturbative part 
of the effective Hamiltonian for heavy-light quark systems, 



(128) 



where 



HQqfreeij = [2(27r)^]^(5^(A;i - k2)S%Pi-p2)Sx,xJx,x, 



x{ {y-l)A + 



—2 -I- „ A 



(129) 



VQqij (y - y', Ki_ - = 2{2TrfS^{ki + pi - k2 - P2)Sx^xJ. 



A2 A4 



-252(T«)(T'^) 



(i^+)2 (y - y')^ 



^ {k± - n'^y _ k]_- K±- k'j_ 

{y-y'Y y{y-y') 



K± ■ k'j_ - (K')i 



y'iy - y') 

0{A- X'^)0{A- B) 
(kx - 02 + {y- y'ft 



0{B - X^)0{B - A) 
(«:^-«l)2-(y-t/')(f -^) 



(130) 



with B 



y-y' 



and the function A has the same form as in 



quarkonium case. Here we have also introduced y = p'l / , k±_ = — yK^, 



but the range of y is now given by < y = < 00. 

The heavy-light quark effective Hamiltonian is mg-independent. This is be- 
cause in heavy-light quark systems the heavy quark kinetic energy can be treated 
as a perturbative correction to Hxq- Obviously the above Hxo has the heavy 
quark spin and flavour symmetry. Compared to the Vq-q, Vq-^ interactions are 
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much more complicated. But it is not difficult to check that the above Vqq con- 
tains a confining interaction. The confining mechanism is the same for QQ and 
Qq systems, as well as for qq systems, as one can show. 

In conclusion, we have obtained the nonperturbative part of a confining QCD 
Hamiltonian for heavy- heavy and heavy-light quark systems. We are now ready 
to solve heavy hadron states on the light-front and to show how the WCT scheme 
works in the present formulation 

5.4 Heavy Hadron Bound States 

As we mentioned the ideas of WCT to nonperturbative QCD is to begin with the 
effective QCD Hamiltonian Hx = Hxo + Hxi- Then using the constituent picture 
to solve nonperturbatively the hadronic bound state equations governed by Hxo 
and to determine the running coupling constant gx- If one could properly choose 
the nonperturbative Hxa such that gx is arbitrarily small, then the corrections 
from Hxi could be computed perturbatively, and we would say that a WCT to 
nonperturbative QCD is realized. Now, I shall discuss such a WCT to heavy 
hadron bound states. 

i). Heavy Hadron Bound State Equation Under WCT. As we have pointed 
out in the first lecture solving eq.(14) from QCD with the entire Fock space is 
impossible. A basic motivation of introducing the WCT scheme is to simplify the 
complexities in solving the above equation. In the present framework, Hlf = 
Hx, where Hx has already decoupled from high energy states. Furthermore, the 
rcscparation Hx = Hxa + Hxi is another crucial step in WCT, where only Hxo 
is assumed to have the nonperturbative contribution to bound states through 
eq.(14), and Hxi is supposed to be a perturbative term which should not be 
considered when we try to solve eq.(14) nonperturbatively. 

The next important step in the WCT scheme is the use of a constituent 
picture. The success of the constitiicnt qiiark model suggests that we may only 
consider the valence quark Fock space in determining the hadronic bound states 
from HxQ- In this picture, quarks and gluons must have constituent masses. This 
constituent picture can naturally be realized on the light-front [Wilson et al. 
(1994)]. However, an essential difference from the phenomenological constituent 
quark model description is that the constituent masses introduced here are A 
dependent. The scale dependence of constituent masses (as well as the effective 
coupling constant) is determined by solving the bound states equation and fit- 
ting the physical quantities with experimental data. But for heavy quark mess, 
this A-dependence can be ignored. Once the constituent picture is introduced, 
we can truncate the general expression of the light-front bound states to only 
including the valence quark Fock space. The higher Fock space contributions can 
be recovered as a perturbative correction through Hxi- Thus, eq.(9) for heavy 
quarkonia can be approximately written as: 



Al A2 
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x<^QQ(t/, K^)bl{h,Xi)dU{k2, A2)|0), (131) 

where the wavefunction (j)^-^{y, kj_) may be mass dependent due to the kinetic 
energy in Hxo [see (117)] but it is spin independent in heavy mass Hmit. Also 
note that the heavy quarkonium states in heavy mass hmit are labelled by the 
residual center mass momentum K^. We may normalize eq.(131) as follows: 

{W{K'+,K'^,X'^)\W{K+,K^,X,)) = 2{27rfK+6%K-K')Sy^x., (132) 

which leads to 



/ 



With the above analysis on the quarkonium states, it is easy to derive the 
corresponding bound state equation. Let Hlf = Hxo of eq.(116), eq.(15) is 
reduced to 



2-t-A. 

rriQ 



2kI + A (2y2 - 2y + 1)J | cj>^-^{y, k^) 



" (..-.i)^!(.-.o-^- ^^"-^^^}^-^^^^'--^-^"^^ 

This is the light-front bound state equation for heavy quarkonia in the WCT 
scheme. 

For the heavy mesons containing one heavy quark, similar consideration leads 

to 

/— 2 —2 + m'^, (A) \ 

(/+ (l-y)/-^_^j<?Q5(y,fc^,Ai,A2) 

2 

= (-l^^'C-zlne) <pQ^{y,k^, Xi,X2) 

^^^Q^(y - y', >^± - '^'±)^Q-,{y', i^'^, Ai, A2), (135) 

where Vq-g is given by eq.(130). Note that the light antiquark here is a brown 
muck, a current light antiquark surround by infinite gluons and qq pairs that 
results in a constituent quark mass rUq which is a function of A. 

ii). A General Analysis of Light- Front Wavef unctions. A numerical compu- 
tation to the bound state equations, eqs.(134) and (135), is actually not too 
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difficult. However, to have a deeper insight about the internal structure of light- 
front bound states, it is better to have an analytic analysis. For this propose, we 
would like to present a general analysis of light-front hadronic wavefunctions. 

The heavy hadronic wavefunctions in the heavy mass limit are rather simple. 
First of all, the heavy quark kinematics have already added some constraints on 
the general form of the light-front wavefunction (j){x,K±). When we introduce 
the residual longitudinal momentum fraction y for heavy quarks, the longitudinal 
momentum fraction dependence in (p is quite different for the heavy-heavy, heavy- 
light and light-light mesons. 

For the light-light mesons, such as pions, rhos, kaons etc., the wavefunction 
4)qq{x, Kj_) must vaulsh at the endpoint a; = or 1. This can be seen from the 

kinetic energy term in eq.(15), where Mg = ^-^^"'^ — for the valence Fock 

space. To ensure that the bound state equation is well defined in the entire range 
of momentum space, \(l)gg(x, kj_)\'^ must fall down to zero in the longitudinal 
direction not slower than 1/x and 1/(1 — a:;) when x ^ and 1, respectively. In 
other words, at least (j)qq{x,Ki_) ~ i/x(l — x) . For heavy-light quark mesons, 
namely the B and D mesons, the wavefunction (f>Qq{y, kj_) is required to vanish 
at 2/ = 0, where y is the residual longitudinal momentum fraction carried by 
the light quark. This is because the kinetic energy in eq.(135) only contains a 
singularity at y = 0. On the other hand, since < y < oo, (f)Qq{y, k±) should also 
vanish when y oo. Hence, a possible simple solution is 4>Qq{y^ '^i) ~ v^^~"^ 
or y^e~"^^. For heavy quarkonia, — oo < y < oo, the normalization forces 
<t>QQ{yi i^i) to vanish as y — > ±oo. Thus, a simple solution may be ({'qqiv^ i^i.) ^ 

On the other hand, the transverse momentum dependence in these light-front 
wavefunctions should be more or less similar. They all vanish at k± — > ±oo. A 
simple form of the k± dependence for these wavefunctions is a Gaussian function: 

The above analysis of light- front wavefunctions is only based on the kinetic 
energy properties of the constituents. Currently, many investigations on the 
hadronic structures use phenomenological light-front wavefunctions. One of such 
wavefunctions that has been widely used in the study of heavy hadron structure 
is the BSW wavefunction [Wirbel et al. (1985)], 



i>BSw{x, K±) = Afy/x{l - x) cxp ( exp 



-^(x-xo)^ 



(136) 



where A/" is a normalization constant, w a parameter of order Aqcd, ^jq = (| — 

2 _ 2 

—2mt^)i and Mh, mi, and m2 are the hadron, quark, and antiquark masses 
respectively. In the heavy mass limit, the BSW wavefunction can be produced 
from our analysis based on the light-front bound state equations. 

Explicitly, for heavy-light quark systems, such as the B and D mesons, one 
can easily find that in the heavy mass limit, mi = niQ ^ Mh, niq « tuq 
so that xo = . Meanwhile, we also have Mhx = Ay. Furthermore, the factor 
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\/x{l — x) can be rewritten by ^ in according to the corresponding bound 
state equation discussed above. Thus, the BSW wavefunction is reduced to 



This agrees with our qualitative analysis given above. Indeed, using such a 
wavefunction we have already computed the universal Isgur-Wise function in 
B ^ D,D* decays [Cheung et al. (1995)]: 

^{vv') = ^, (138) 
V ■ v' 

and from which we obtained the slope of ^{v ■ v') at the zero-recoil point, = 
— ^'(1) = 1, in excellent agreement with the recent CLCO result [Patterson 
(1995)] of /a^ = 1.01 ± 0.15 ± 0.09. 

For heavy quarkonia, such as the hh and cc states, mi = m2 = mq which 
leads to xo = 1/2 in eq.(136). Also note that Mh{x — 1/2) = Ay, and the factor 
\Jx{l — x) must be totally dropped as we have discussed form the quarkonium 
bound state equation. Then the BSW wavefunction for quarkonia is reduced to 




<^QQ(y,'ti) =^ exp (^-^j exp ( ~—y' ) , (139) 

which is a form as we expected from the qualitative analysis. Here we have not 
taken the limit of mg 00 for heavy quarkonia. Thus a possible toq dependence 
in wavefunction may be hidden in the parameter w. 

Using the variational approach with the above trial wave function, we can 
analytically solve the bound state equation Eq.(134). Under the consideration 
of SRG invariant for the binding energy A, we find that the effective coupling 
constant [Zhang (1996)]: 



(140) 




a-|-61n^ 

A 



where the coefficients a and h are determined by minimizing the binding energy 
from (134) with (139). The coefficient b is almost a constant (with a weak de- 
pendence on vriQ but independence on A and A), while a depends on both A and 
toq, and also slightly on A. For A > 0.6 GeV, the A-dependence in the parameter 
a is negligible. It is known that A is of the same order as Aqcd which is about 
100 ~ 400 MeV. The charmed and bottom quark masses used here are mc = 1.4 
GeV and mb = 4.8 GeV. From the particle data [Particle Data Group (1994)], 
the binding energy A should also be less than 400 MeV. For a qualitative con- 
sideration, we may take A = 0.2 GeV and A = 0.4 GeV for charmonium. Then 
we find that b = 1.15, and a = —0.25 and 1.1 respectively. 
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(141) 



which is imich smaller than that extrapolated from the canonical running cou- 
pling constant in the naive perturbative QCD calculation. In order to see how 
this weak coupling constant varies with the scale A, we take A = 200 MeV 
and vary the value of A around 1 GeV. Wc find that the coupling constant is de- 
creased very faster with increasing A. In other words, with a suitable choice of the 
hadronic mass scale A in SRG, we can make the effective coupling constant a\ 
in H\ arbitrarily small. Then the WCT of nonperturbative QCD can be achieved 
in terms of Hx such that the corrections from Hxj can be truly computed pertur- 
batively. This provides the first realization of the WCT to nonperturbative QCD 
dynamics on the light-front. 

5.5 Perspectives 

The applications of the present theory to heavy quarkonium spectroscopy and 
various heavy quarkonium annihilation and production processes can be simply 
achieved by numerically solving the bound state equations (134). and by fur- 
ther including the l/mq corrections (which naturally leads to the spin splitting 
interactions). The extension of the computations to heavy-light quark systems 
is straightforward. The extension of the present work to light-light hadrons re- 
quires the understanding of chiral symmetry breaking in QCD which is a new 
challenge to nonperturbative QCD on the light-front. Nevertheless, we have pro- 
vided a detailed analysis to the weak-coupling treatment of nonperturbative 
QCD proposed recently [Wilson et al. (1994)]. I believe that LFQCD opens a 
new research direction in the attempt of solving the most difficult problem in 
field theory, that is, the problem of the rclativistic composite particle bound 
states governed by the nonperturbative dynamics of QCD. 
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